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Abstract 

The screening of magnetic Z(A)-monopoles and the associated screening length 
in SU(N) gauge theories are analyzed theoretically, and computed numerically 
in the 3d SU (2) theory. The nature of the screening excitations as well as their 
mass have so far remained inconclusive in the literature. Here we show that the 
screening mass is identical to the lowest J^'-'' = 0^^ excitation of the 4d Yang- 
Mills Hamiltonian with one compact direction with period 1/T, the subscript 
R referring to parity in this direction. We extend the continuum formulation to 
the lattice, and determine the transfer matrix governing the decay of the spatial 
monopole correlator at any finite lattice spacing. Our numerical results for 
SU{2) for the screening mass in the dimensionally reduced (high temperature) 
theory are compatible with the 0"*""^ glueball mass in 3d SU{2). 



1 Introduction 



Ever since the recognition by 't Hooft and Mandelstam pP that confinement in QCD 
might be a dual version of superconductivity, people have sought for a quantitative 
understanding. The first step in that direction was the proof that electric confine- 
ment implied the screening of colour magnetic fields j2] . That is of course a necessary 
condition: the would-be magnetic Cooper pairs are expected to screen colour mag- 
netic fields. The proof in ref. |2I also pioneers the method that we will use below 
to extract the magnetic screening mass for any temperature, namely by introduc- 
ing heavy magnetic sources and measuring their correlation, 't Hooft used periodic 
boundary conditions and introduced sourceless fiuxes winding around the periodic 
space directions. Electric confinement in this terminology means the electric fiuxes 
carry an energy F^. ~ La. Here cr is the string tension per unit length and L the size 
of the box in the direction of the fiux. Now one makes the mild assumption that at 
zero temperature the energy of any array of electric and magnetic fiuxes is additive, 
i.e. splits into a part Fg due to electric fiuxes and in a part Fm due to magnetic 
fiuxes. It then follows that F^ ~ crL exp(— crL^^), i.e. decays exponentially with the 
cross section orthogonal to the fiux, the proportionality factor being the energy ^/a 
per unit length. 

Instead of studying periodic fiux loops, one may gain more information by looking 
at a monopole-antimonopole pair along the z-direction and by varying its separation 
r, in general at some finite temperature T. This yields a free energy Fm{r,T) as 
a function of monopole separation, which is typically parametrized by a Yukawa 
behaviour. 



Following the lattice formulation of the problem [3 H], there have been many at- 
tempts to compute the screening mass rriM by numerical simulations |3]-[TT], and 
to interpret the corresponding excitations as effectively massive 'constituent' gluons, 
Debye-screened gluons or glueballs. Since an accurate numerical computation of the 
screening mass tum is very difficult, the numerical values are rather scattered and 
these studies have remained inconclusive. 

In this paper we show that the screening masses rriM correspond to spin zero mass 
levels with PC = 1 of a spatial Hamiltonian Ht- This Hamiltonian propagates gauge 
invariant 3d Yang-Mills states (i.e. glueballs) along the ^-direction. The subscript T 
means this Hamiltonian describes 3d Yang-Mills fields in a box with the z direction 
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periodic with period ^ and the x and y directions infinite in extent. So the rotation 
group in the x — y-plane is SO {2). Let parity P be refiection in the two-dimensional 
infinite space, and R be the refiection in the periodic direction. Together with charge 
conjugation C the corresponding quantum numbers are conserved by Hrp- Until now 
a quantum number assignment for the screening states was not mentioned in the 

literature, as far as we know ^. Here, we find that the assignments for the magnetic 

P c 

' are J = and R = CP = 1. 

Our underlying reasoning is similar to that applied to electric screening by Arnold 
and Yaffe some time ago ^21- The operator creating the heavy monopole is a local 
operator, in the same sense as the operator creating a heavy quark. The locality 
of the operator is somewhat masked in the path integral representation, which may 
explain why it remained unnoticed. For this reason we will be fairly explicit. We 
discuss the lattice version of the operator in detail and establish that the relation 
between the screening mass ttim and the mass of the lowest 0^^ state is true for 
any finite lattice size T~^. Two limiting cases are worth mentioning. At T = the 
screening mass niM corresponds to a 4d scalar glueball with R = P because of full 3d 
rotation symmetry. Since CP = i? = 1, we thus get the O"*""^ glueball mass. At large 
temperatures, modes of order ~ T can be integrated out and screening masses are 
well described by the dimensionally reduced theory (see ^H! and references therein). 
The latter corresponds to a Yang-Mills theory for the colour magnetic fields Ai ~ g'^T 
coupled to an adjoint scalar Aq ~ gT in 3d. At T » mu, the Debye screening mass, 
the electric fields Aq ~ gT may be integrated out as well, and our screening mass 
corresponds to a 0^"-^"^ state of the 2d Yang-Mills Hamiltonian, i.e. a glueball of the 
3d Yang-Mills theory. 

This connection is of great numerical value, since simulating the monopole pair 
is numerically quite involved |H1 E]. Exploiting the connection, one might instead use 
a simple local source with the same quantum numbers as the heavy monopole and 
measure its correlation by well-established ^H], numerical methods in a Yang Mills 
system at zero temperature and with one periodic direction. 

The layout of the paper is as follows. In section |H monopole sources are discussed. 
Their correlation and the connection to the spatial Hamiltonian is the subject of 
section IHl and section |31 In section we will discuss the lattice operator that excites 
the state with the screening mass, and the relation of its correlator with the object 
we simulate: the twisted action. We identify the transfer matrix whose lowest scalar 

^Some of it is discussed by one of us in the 2003 Zakopane lecture notes |14| . 
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glueball mass coincides with the screening mass from our lattice correlator. Finally, in 
sectionElwe will put this connection to the test by simulating at very high temperature 
the thermally reduced version of SU (2) gauge theory. The restriction at very high T 
permits us to simplify the monopole correlator to one Euclidean time slice, i.e. to a 
two-point funtion. Our results compare very well with the known masses ^HI in the 
corresponding 2d reduced Hamiltonian. Readers who are interested in the numerical 
results only may want to skip the first sections. 

Before delving into the technical part of the paper a word on motivation. The 
reader might wonder whether the study of magnetic screening is an entirely academic 
exercise. We do not believe so. Understanding magnetic screening quantitatively is a 
necessary step in the understanding of the magnetic activity that causes it. And this 
magnetic activity is at the core of confinement 

Apart from this matter of principle there is a computational reason already partly 
mentioned before. By rotating the time direction into a space direction one gets a 
spacelike 't Hooft loop, or "electric flux loop" , with an area law behaviour. The loop 
tension does not correspond to a level in the fictitious Hamiltonian Ht, although the 
simulation is identical in complexity [SI E]- Thus a calibration of these simulation 
methods is more than welcome and that is precisely provided by the connection of 
the magnetic mass with a mass level accessible through usual methods . 

2 Creating a heavy monopole 

In this section the analysis of the creation operator of a Dirac monopole is resurrected. 
It may well have been done way back in the seventies, but we have been unable to 
find a reference. To a large extent it is related to the vortex operator in the seminal 
work of 't Hooft in 1978 [2], and used extensively in the work by Kovner J7]. We 
apologize for the pomp and circumstance that will go with this resurrection, but we 
believe it is useful to insist on the basic physics involved in order to digest the sections 
that follow. 

We consider an SU{N) gauge theory in four dimensions, and matter fields that 
have A^-ality zero. So they are neutral under the centergroup Z{N) of SU{N). Con- 
sider in such a theory the correlator of two heavy "electric" colour charges Qk with 
non-zero N-ality k. Qk is the creation operator of a particle field in a representation 
of SU(N) with N-ality k and with an infinite mass. A centergroup transformation 
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exp(i^) transforms the source as follows : 

27r 

Qk exp{ik—)Qk. (2) 

The correlator Tr (^^(O, 0, r) exp(— if/T)(5fc(0) in path integral language corresponds 
to the well known correlator of Polyakov loops, 



P(A^(f)) = Pexp ii Ar\ 



X, T)dT , (3) 



up to obvious normalization factors: 

Trgl(0,0,r) exp{-H/T)Qk{0) = J DA P^A, {0,0, r))P{Ar (6)) exp (^-^S{A)j . 

^ (4) 

At- is in the same representation as Q^. The exponential decay of such a correlator, 
the static potential, is used as an order parameter: at low temperature the potential 
rises linearly, at high temperature it is screened like in eq. (jT]). 

Here we want to formulate its magnetic analogue. In this section we first discuss 
the heavy monopole source, before using it to construct the magnetic correlator in 
the next section. 

2.1 Dirac monopole in electrodynamics 

The monopole source in SU{N) will be like the Dirac monopole in electrodynamics, 
which we will discuss first. The simplest picture of a magnetic charge at the origin is 
one with a radial magnetic field 

Bmon = (5) 

yielding a total charge m from our monopole after integrating over any surface around 
the monopole. It is a solution of the Maxwell equation V.Bmon = m6{x). 

Unfortunately the source term on the r.h.s. excludes the use of vector potentials 
because of the Bianchi identity, which says that a magnetic field expressed in terms of 
a gauge potential has no divergence. Dirac solved that problem by putting the source 
term at infinity. He then put a string of magnetic dipoles from that source to x = 0, 
guiding the flux through this thin string from the monopole at infinity to x = 0. Let 
this magnetic "return flux" come in along the z-axis T^. Then it is given by: 

Bs = VxS= -m6{x)6{y)e{z%. (6) 
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Adding to the monopole field strength the return flux field strength, 

Bmon + Bg = Bsol, (7) 

gives a configuration with divergence zero everywhere. It is the field strength of a 
solenoid along the positive axis. By taking the surface integral it becomes explicit 
that the solenoid has no magnetic charge, because the integrated fiux of the monopole 
field is neutralized by incoming fiux of the string. So there must be a potential Agoi 
with V X Agoi = Bsol- Its precise form is not relevant for what follows. 

Hence the string represents the source of magnetic charge just as the delta function 
represents the source of electric charge. The end point of the string is where the 
magnetic charge resides. If we are interested in having two monopoles of opposite 
sign then the string does not come in from infinity, but from the location of the second 
monopole with the opposite charge, and runs, as before , to the monopole eq. © at 
the origin. 

Of course the string should be invisible by scattering with a quantum mechanical 
particle with charge e. Then the string with its vector potential written as Ag causes 
a phase difference 

ex.p{ie (j) dl.Ag). (8) 

Now Stokes' theorem § dl.Ag = J dS.Bg = m tells us, that for the phase difference to 
vanish, the Dirac condition |flB- em = 2Tm must apply. 

Finally we want to write down a monopole creation operator in QED. As we just 
saw, this amounts to creating an invisible string, obeying the Dirac condition. The 
string has magnetic field zero everywhere except on the string itself, running along 
the positive z-axis. Thus the vector potential Ag of the string is a pure gauge Vuj 
everywhere, except on the positive z-axis. There the gauge transformation is singular 
in such a way as to recover the Dirac condition eq.®. For any closed circuit around 
the string the gauge transform uj has a discontinuity Auj with: 

e j dl.Ag = Auj = k27i. (9) 

So taking u = k(j), the azimuthal angle around the string, we see that going around 
the string once, from = to = 27r, the gauge transformation is forced to go k 
times around the circle that constitutes the gauge group U{1). Such transformations 
have a non-trivial homotopy Hi. They are written as uJk{xo]x), xq being the point 
where the string starts and runs along the positive z-axis. 
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Thus, to create a heavy monopole source in QED we have to effect a singular gauge 
transformation, generated by the Gauss operator G{x) = V.E — ep, with charge 
density ep. A physical state in the Hilbert space is by definition invariant under 
regular transforms, in particular with trivial homotopy. So the monopole operator 
Mi(xo) with charge = 1 in units of ^ equals: 

Mi{xo) = exp J dx uJi{xo; x)G{x)^ . (10) 

2.2 Monopole source in SU(N) gauge theory 

Let us now turn to the generalization of the Dirac monopole to the SU(N) case. We 
take gluodynamics, possibly with matter fields in the adjoint representation. The 
strength of the SU(N) Dirac monopole is again given by the discontinuity of the 
gauge transform encircling the string, the discontinuity now given by a centergroup 
element. Because we suppose absence of charged Z(N) fields, the gauge group is 
really SU{N)/Z{N). Encircling the string gives no singularity in this group, but 
corresponds to a non-contractible path, just like for U{1). 

Having identified the subgroup Z{N) we now proceed to explicitly construct the 
operator that creates the Z{N) monopole. We need to describe the string emanating 
from the monopole at the point 0, say in the positive 2;-direction. For that we need 
a gauge transformation Q{x) that has a discontinuity in the centergroup, when going 
around the string. To this end we define a function in the Lie-algebra of SU{N): 
u{x) = (j){x)^Yk. X is any point on the string and 0(x) is the azimuthal angle when 
following a closed curve around x in the x — y plane. Yk is a special traceless diagonal 
matrix in the Lie-algebra, which is the "k-hvpercharge" . = diag(/c, k, k, k — 
N, k — N, ...k — N), such that when exponentiated it gives a Z(N) group element: 

exp (i'^Yk^ = exp (ik'^^ Inxn- (H) 



N J " V 
This means the diagonal SU{N) gauge transformation 



fifc = exp (^iu;(x)^Yfe^ (12) 

picks up a discontinuity exp{ik^), when encircling the string. 

The field strength of the string for the SU(N) Dirac monopole then reads BgYk 
in the fundamental representation, with Bg given by eq.©. There is an important 
question left. Under regular gauge transformations this field strength transforms as 
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Since a regular gauge transformation does not change the discontinuity, one would 
expect the total flux $ of the monopole not to change. That means that the only 
gauge invariant flux one can have is exp{imYk) with m, the strength of the unit 
magnetic flux through the string, equal to 

This motivates the definition of the Z(N) Dirac monopole operator with charge 
exp(i/c^) as the operator representation of the gauge transformation in eq.lfT^: 



The dot means integration over the space coordinates and a trace over colour. 

The string will not affect physical states, built from Z{N) invariant matter. Phys- 
ical states will only contain Wilson loops in Z{N) neutral representations, and these 
loops will not respond to the string piercing them. Therefore the location of the string 
is immaterial under these circumstances, only its endpoint matters. This means the 
monopole source is a scalar under rotations. Under P and C the orientation of the 
string reverses. The gauge invariant magnetic flux changes therefore from exp(z^) 
to exp(— i^). The rotations and parity discussed here are those of the 3d-space in 
which we live. In section |3] we will discuss the parity P and R parity in the fictitious 
3d space with one direction compactified, and construct an operator that creates the 
corresponding state. 

At this point it is clear that the correlation of two of our static magnetic scalar 
sources is analogous to that of two static electric sources as in eq. (jl)): we are looking 
at the magnetic Coulomb force and its screening. Note that this screening mass is 
different from the pole mass of the magnetic gluon propagator considered in other 
work [20!. 



3 Correlator of two Z(N) monopoles at finite T 



Now that we have found the local operator that creates the monopole, we want to 
know the force law between two of them. To formulate the force law, we start with 
a monopole at 0, an anti-monopole at (0,0, r) and form their correlation. This is 
the operator Mk{0, r) given by a gauge transformation having the now familiar Z{N) 
singularity on a string going from to (0, 0, r). Running it in imaginary time gives: 



for the vacuum to vacuum amplitude to create and annihilate our monopole pair. 
F(r, r) for r — > oo is the energy associated with the pair. 




(13) 



C(r) =< 0|Mfc(0,r)exp(-iJr)Mfc^(0,r)|0 >= exp(-rF^(r, r)) 



(14) 
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Figure 1: Monopole antimonopole pair induced by twisting the (x,y) plaquettes 
pierced by the Dirac string in the ^-direction. 

Every Wilson loop with N-ality / encircling the string will pick up a phase exp{ikl^). 
I is the number of fundamental representations the loop is made from. If the Hamil- 
tonian contains a regularized version of the magnetic field strength density, as on the 
lattice, then: 

Bl ^ W^y + Wl- (15) 

This is a somewhat hybrid method, half continuum, half lattice, and we will give 
a full lattice version in section El The Wilson loop W^y is taken in the fundamental 
representation / = 1. Then it is clear that for a fixed time slice Ar one can commute 
the correlator through the Hamiltonian factor exp(— ifAr) and produce a twisted 
Hamiltonian Hk{r), with all loops encircling the string replaced by |21j : 

W.y + Wly ^ exp {ik^^ W^y + exp (-^fc^) Wly (16) 

This is shown in fig. (^. Repeating this for all time slices, and going to the path 
integral representation of eq. ()14|1 will give us the usual path integral in which the 
action S is replaced by the twisted action S'(fc)(r, r)PT|: 

exp(-rF^(r,r)) ~ j Z^Aexp |^-l5(fc)(r, r)^ (17) 

The xy plaquettes on the sheet traced out by the string are all twisted like in eq. ()l(j|l . 
This creates a (2;, r) temporal 't Hooft loop. 

The extension to finite temperature T is now simply: 

exp (^-lF^(r,r = 1/T)) ~ j DA exp (^-l5(fc)(r, r = 1/T)j (18) 

where the sheet now wraps around the full periodic time direction. In the following we 
shall drop the argument 1/T in the action. The free energy Fm{r,T = 1/T) depends 
on T and is, as discussed above eq. ((H), screened for all temperatures T. 
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4 Space-periodic Hamiltonian and its symmetries 



The correlation in its path integral form eq. ()18|) can be read as a propagator in the z- 
direction propagated by a Hamiltonian Ht- The "time" direction for this Hamiltonian 
is the 2;-direction. Its space directions are the x, y and periodic r direction. It 
correlates a "magnetic Wilson line" wrapping around the r-direction at 2; = 
with another one at z = v. 

CiU){r) = j DAeM-^S^k){r)) 

~ Tr exp(-iJr(L-r))\4exp(-/7Tr)Vfc^ (19) 

where L ^ 00 and the trace is over physical states. The "magnetic Wilson line" 
is the operator 

Vk = exp (^t fj^ dTE.Duj{x, y) , (20) 

which is the analogue of the "electric" Wilson line P{Ar). The tilde indicates that our 
canonical operators are now At,x,?; and ET-^x,y The gauge function equals u!{x,y) = 
arctan(?//x). Hence this operator creates for every value of r a vortex of strength 
exp(z/c^) in the origin of the x, y plane. 

The Hamiltonian Ht describes the Yang-Mills field in two infinite and one periodic 
dimensions and propagates it in the ^-direction. For the 2;-slice Az between z = 
and z = 1 the line operator Vk will twist the Bz operator in this Hamiltonian over the 
full period 1/T in r. This can be seen in the same fashion as in the previous section 
by commuting the Vk operator through the factor exp{—HTAz). Repeat this for the 
second slice and so on. We reproduce the same sheet of twisted x — y plaquettes as 
in the previous section, hence the identity eq.(fT^. 

Inserting a complete set of states tells us that the lowest mass in the set of 
intermediate states will be the screening mass^: 

<0\Vkexp(-HTr)VhO> = I < OlV^lO > 1^ + I < 0\Vk\mM > |2 ^^P(~"^^^) ^ 

r 

= exp(--Fo)(l + Cfc + ....). (21) 

I r 

If we want to create states with all momentum components zero, we have to inte- 
grate the sources over x and y directions as well. This however would not correspond 
^Normalization is such that the energy of the vacuum state is zero. 
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to what we are doing on the lattice. Our source is fixed in x and y. Thus, in identi- 
fying from eq. ()2H1 the lowest mass, we have to integrate over momenta in the large 
directions x and y. So that is why in that formula we have a Yukawa potential: 

I < om,n,, > r ^^P'-'"'"-' . (22) 

r 

resulting from the integral / ^^^^^ exp(— £'(p)r). Neglecting momentum dependence 
in the matrix element is allowed at large enough r. Now we have to identify the 
quantum numbers of Vk- 

4.1 Quantum numbers 

Consider the symmetry group ^ of Ht'- S0{2)xPxCxR. Apart from the 2d rotation 
group there is 2d parity P, which fiips the sign of the y-axis, and Ay —>■ —Ay. There 
is charge conjugation C with Ar,x,y — ^ —Ar,x,y, and R parity changing t ^ ^ — t, 

Clearly Vk couples to scalars under SO (2), since it is running parallel to the SO (2) 
rotation axis. Under P and C it transforms into V^, but R parity leaves it invariant. 
So our magnetic screening mass is a 0^'-' state with PC = R = 1. For the symmetric 
combination Vk + we have P = C = 1, for the antisymmetric combination P = 
C = —1. Periodicity in the number of colours and charge conjugation tell us that Vk 
and V/v-fc have the same effect. 

However we can fix the assignment completely by reducing the temperature to 
T = 0. Then we restore the full rotational symmetry. Since R and P refiections are 
related by the generator rotating the y-direction into the r direction, a spin zero 
state like our magnetic mass state at T = must have R = P. That excludes the 
anti-symmetric combination mentioned before and only 0''"^ is possible. 

An amusing question comes up in connection with the correlator of Im V^. If it is 
to be different from the correlator of ReV^, the correlator of Vk should be non-zero. 
The question then arises how to implement the twisted path integral version of such 
a correlator. Clearly one needs the return flux from a heavy monopole at inflnity 
with strength —2k (or N — 2k). To simulate such a twist conflguration is in principle 



■^In fact this is not all of the symmetry group of Ht- There is the Z{N) group generated by 
gauge transformations periodic modulo Z{N) in the r direction. They have as order parameter the 
Wilson hne in the r direction, P{Ar). Below Tc we have < 0|P(At-)|0 >= whereas above Tc Z{N) 
is spontaneously broken: < 0|P(ylT-)|0 >^ 0. Our operator Vk does not transform under Z{N), nor 
any local state. Only P{At) does. 
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doable, but hard in practice ^. But in practice one would -and could, as we just learnt- 
excite the same states by some local operator with the same quantum numbers as 

It is instructive to compare with the electric Coulomb force, given by the Wilson 
line P{Ar{x)) = TrPexp(i Jq^^ drAriT, x)). C and R change it into P{Ar)\ while 
P leaves it invariant. So it excites 0^*" states with CR = P = 1^ C = R = 1 
for the symmetric, C = R = —1 for the antisymmetric combination, and radial 
excitations. Like in the magnetic case we can narrow down the quantum numbers, 
as argued in ref. ^2)- For small enough coupling we can use one loop perturbation 
theory, and the self-energy of Aq is clearly odd in R. Note a curiosity in the SU{2) 
theory. There the imaginary part of P{Ar) is zero, so it cannot excite states with 
R = —1. Nevertheless states with R = —1 do exist in the spectrum of the periodic 
Hamiltonian H^l Their mass has been computed including the first non-perturbative 
contribution |22[ IT8]. However, for all SU{N) theories the correlator of lmP{Ar) is 
exponentially suppressed below (like exp(— a being the string tension). 

4.2 Dimensionally reduced monopole correlator 

What happens under dimensional reduction to our monopole correlator? Obviously 
the one time slice version of our twisted action survives, and all one has to do in a 
3d simulation is to use the correlator along one string stretching from to (0,0,r). 
The free energy F{r) is then again a Yukawa potential (temperature dependence is 
absorbed into the parameters) 

F^{r) = Fo exp(-mMr). (23) 

r 

The answer one gets in the continuum limit is expressed in terms of the 3d coupling 
gM = g^T: 

^ = 4 (24) 

9m 

The dimensionless number d gives the result for asymptotically high T. According 
to the identification of the magnetic screening mass in the previous sections it should 
equal the ratio of the 0^^ mass over the dimensionful coupling g^j in 3d Yang-Mills. 
This is indeed what we find by numerical simulation and is discussed in section IHl 



*We thank Christian Hoelbhng for discussions on this point. 



11 



A 



A / 


A 


A 


i 


■ ■— 1 

/ 1 / 


A 


A 
















I i. 
















V 




V 




A 




A 




A 




A 





Figure 2: Change of the string location brought about by the redefinition of the gauge 
variable on the link denoted by a bolder line. 

5 Lattice formulation 

In this section we give the lattice formulation of the identification of the screening 
mass with the mass of a state of the periodic Hamiltonian. Indeed, as emphasised in 
sectional a consistent formulation requires a regulator for the magnetic field strength 
in the Hamiltonian, and a natural regulator is given by the lattice. In particular 
there is a subtlety with the choice of the transfer operator on the lattice. We will 
need a choice with the same spectrum as the conventional one, which at the same 
time produces the twisted lattice action as defined in sectional 

For simplicity and numerical feasibility, we work in the dimensionally reduced 
formulation. At the end of this section the extension to the 3+1 dimensional case turns 
out to be quite straightforward. We furthermore restrict our theoretical discussion 
to the simpler case where the spatial [x — |/)-extent of the system is infinite. Our 
simulated volumes are chosen large enough so that finite volume effects are negligible, 
as we will argue below. 

Let M be the set of plaquettes being pierced by the Dirac string. They pick up 
a Z(N) factor (^^ = exp(ifc^). This is shown in Fig. 2J Then the lattice action for a 
Euclidean pure gauge theory in the presence of a monopole pair separated by r is 

= 4 E ReTrf/(p) + ^ 5: ReaTrf/(p), (25) 

which is obtained from the standard Wilson action S by simply multiplying all plaque- 
ttes j9 G M by the Z(N) factor C,^. Correspondingly, this defines a partition function 
Z{k){r) as a path integral evaluated with S(k){r). We are interested in the behavior 
under growing r, and thus consider the correlation function 

q,)(r) = = (^e^^..M^<^^-^^^^^M'^ . (26) 
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Before we carry on we note that the location of the string in the path integral is 
not important. This is illustrated in fig. One takes a link variable U{1) on any 
of the pierced plaquettes, and uses the invariance of the measure to change it into 
U{1) = CkU{l). This has the effect shown in the figure: the original twisted plaquette 
loses the twist Cfc when expressed in f/(/), but the other plaquettes bordering our 
link will acquire the twist. A deformation in the string results, but the path integral 
remains the same. Only the end points of the string are fixed as the reader can easily 
verify. This suggests that also on the lattice there is a local operator Vl(xo) at the 
end points, creating and annihilating the vortex. 

For notational convenience we drop from now on the reference to the vortex center 
Xo, and to the subscript L. We will write simply Vk for the operator Vl- In continuum 
language this operator was a singular gauge transformation, and our task is now to 
find its lattice equivalent. More precisely, we want a lattice version of the Hamiltonian 
formula (|T9|l . valid in the continuum: 

C(fc)(r) ~ Tr e^v{-HT^)Vkew{-HTr)Vle^^{-HT{^-r)). (27) 

In a lattice formulation we wish to express the r.h.s. in terms of some transfer matrix 
T' and the operator Vk- In contrast to the standard formulation jSlIll; we choose the 
transfer matrix to propagate states in the z-direction. That is, with the size Lz in 
the 2-direction and distance r in units of the lattice spacing, the correlator can be 
written as 

C(fc)(r) ~ TrT'^r,r''^r;T'(^-'-). (28) 

In this form one has established exponential decay with the spectrum of T' in the 
quantum number sector of Vk-, and the screening excitations are easily identified. 

5.1 The lattice vortex operator 

The vortex operator in eq. (|28|) is given by its action on a state \U{1) >, where inside 
the ket the whole collection of SU (N) link variables in a fixed time plane (with x — y 
coordinates) is given. The operator valued matrix U_ is diagonal in this basis: 

U\U>=U\U>. (29) 

The action of the vortex operator on such a ket is simple. Draw from the vortex 
center in the middle of a plaquette, xq = (xq + a/2,yQ + a/2), a line to infinity, say 
parallel to the x-axis. This line will cross a set of links in the positive ^/-direction. 
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The vortex operator multiplies all those links with the Z(N) phase Cfc = exp(i/c^). 
They are the "twisted" links. The rest of the links stays unaffected. We write in a 
shorthand notation: 

Vl\U >= ICfef/ > . (30) 

As we said explicitly before, the x — y plane is taken to be infinite in extension. This 
avoids problems with the endpoint of the line that cuts the twisted links. As a result, 
the plaquette containing the vortex at its center gets twisted, while all others remain 
untwisted. Our definition tells us that Vk is unitary, just like the continuum operator. 
This can also be seen from the product {U'\V^Vk\U) = {zU'\zU) = 6{CkU\C,kU) = 
6{U', U), since the delta function on the group is invariant under any rotation. 

An explicit representation of the vortex operator can be given in terms of the 
canonical variables conjugate to the U_, the electric field operators defined by their 
commutation relations [£"{1), E\l)] = if^^E^l), [E^{l),U{l)\ = -T^U{1) With 
E°'{l)T°' = E{1), the vortex operator can then be written as 

Vk = Uiexp(t^TTE{l)Yky (31) 

The product is over all links that are crossed by the line emanating from the vortex 
center xq. 

Now all Wilson loops in the fixed 2;-plane that circumnavigate the vortex center 
will pick up the Z(N) phase. Hence, transforms a physical state into another 
physical state. A physical state \ip > is a. linear superposition of Wilson loop states 
Y^l^'lIL >. Applying to such a state will change by a Z(N) phase those loops 
L that circumnavigate the vortex center, creating thereby a new physical state. It 
is then clear that the effect of on a physical state does not depend on the way 
we defined the line of twisted links. This is another way of noting that a gauge 
transformation in one of the end points of a twisted link results in a deformation of 
the set of twisted links. 

We thus expect the vortex operator to be a scalar under the remnant of the 
rotation group admitted by the lattice. Under parity (in d=2: y —>■ —y, x ^ x) and 
charge conjugation we have Vk —>■ Vj. Hence in the continuum limit we expect two 
sets of excitations: a set 0++ and a set of . 

5.2 The transfer matrix 

We define the transfer matrix T = exp(— aif^) to propagate quantum mechanical 
states by one lattice spacing along the 2;-direction, i.e. z plays the role of time in our 
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Hamiltonian treatment. Starting from the conventional definition of T [53], we also 
introduce a more unconventional T' which is needed in order to write our correlator 
in the form eq. ()28p. The difference between the two transfer matrices is in their 
decomposition into kinetic {K) and potential {V) parts as 

T = Ty'^TxTy'^ (32) 

r = t]I^TvT]I^. (33) 

Before defining the factors, we note that T and T' have the same spectrum, and hence 
they yield identical partition functions. This is inferred from the fact that the trace 
of any power of T and T' is the same, by using cyclicity of the trace: 

Tr = Tr T'^ = Tr (TvTk)^ for any N. (34) 

The different factors are constructed from the potential and kinetic terms of the 
action, respectively. The potential factor Ty is a multiplication operator and consists 
of the exponent of the sum of all spatial (x — y) plaquettes pij: 

Ty = exp (j^ Re Tr L/(p,,) j . (35) 

This operator is diagonal in the [/-basis, \U). The kinetic operator is an integral 
operator and defined by its matrix elements (r = {x,y),Vi unit vector in i-direction, 
i = x,y): 



< U'\Tk\U > 



\ r;i=x,y j 

= J DU,exp(j^J2^eTTU{poi)^ (36) 

On the r.h.s. of this definition appear all the "time"- like plaquettes poi ( i-e. with 
one link in the z-direction). In temporal (viz. axial) gauge {11^ = 1 on all z-links), 
this definition reduces to: 



< U'\Tok\U >= exp Y ReTrf/,(f)f/;^(f)) ] . (37) 

r;i=x,y 



N 



The operators in axial gauge and the gauge invariant ones are connected by a projec- 
tion operator P, 

Tk = TokP = PToK, P = j ^9 R[gl (38) 
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where the unitary operator R[g\ induces a gauge transformation with gauge function 
g, R[g]\U>=\W>. 

Ty is a multiphcation operator in the U-basis, so its square root is simply defined 

as: 

= """P (^EReTrI/(p,,) J . (39) 

\ Pi] I 

On the other hand, from our definition of Tk and Tqk it is not immediately obvious 
how to define their square root, since they are non-diagonal in the [/-basis. A pos- 
sibility is to define T^j^ by its matrix elements in terms of the character expansion. 
The character Xr{U) is the trace of the (i^^- dimensional representation matrix D^{U) 
in the given representation R, Xr{U) = TrD^{U). The character expansion of the 
kinetic factor is 

(U'ITokIU) = Y.dRXR{UU'^)lR{l3). (40) 

R 

Assuming that the square root is also a function of ReTr {U {li)U'\li)) , its character 
expansion is 

{U'\{TokY"\U) = Y.dRXR{UU'^)JR{(3), (41) 

R 

with as yet unknown coefficients Jr. To find these we note that the square root of 
the kinetic factor is required to satisfy 

{T,kY'\T,kY''=T,k. (42) 

In terms of the matrix elements this means 

' dU"{U'\{T,KY''\U"){U"\{T,KY''\U) = mTm\U). (43) 



Using eq. (|¥T|) and the orthonormality of the characters ^, 

/ dU"xR{UU")xR'{{U'U"Y) = 5a,R,^XRiUU'^), (45) 
one finds that eq. PHj) reduces to: 

{U'\Tok\U) = J2dRXRiUU'^)JU(3). (46) 

R 



^This formula follows from the orthogonality of the representations: 

/ dUD'^iU)m,nD'^'{U^)n',ra' = 6 R,R, dm,m' Sn,n' ■ 

J an 



(44) 
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Comparing coefficients with eq. ()40p one determines 

JniPf = W). (47) 

Now we define 

T'J' = TiL'P, (48) 
and observing that T^I^P = PTqj^, one easily verifies that TjJ'^TjJ'^ = Tk as desired. 

5.3 The twisted transfer matrix 

Next, we consider the action of the vortex operator on the transfer matrix. Acting 
on the kinetic operator in axial gauge, its matrix elements are 

< U'\VkToKVi\U >=< CkU'\ToK\CkU >, (49) 

according to the definition of Vk- Using eq. (jHZj) we get 

< CkU'\ToK\CkU >= exp^{Y^ReTTiCkU'iQyCkUilt) + J2^eTTU'il'J^Uilnt)}. 

h Int 

(50) 

Only in the twisted links (It) there is the centergroup factor, not in the untwisted 
links Int- Of course all the centergroup elements drop out of the r.h.s. of eq. ()5U|) . 
so VkToxV^ = Tqk ■ This is true for any lattice action defined in terms of a sum 
of traces of irreducible representations. Furthermore, the vortex operator commutes 
with the projection operator P, since the latter does not contain any link variables. 
Hence we also have 

VkTKV^ = PVkToKV^ = PToK = Tk. (51) 

On the other hand, the potential factor of the transfer matrix contains the pla- 
quette which encircles the vortex center, and this plaquette picks up a twist. We 
therefore obtain the twisted potential operator 

Tf = VkTyV^ = ( 1^ ^ ReTr f/(p,,) I exp ReCfcTr t/(p,,(fo))) • (52) 

With these results the twisted transfer matrices take the form 

T'W = VkT'V^ = T]I^T^^T]I\ (53) 
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We are now able to express our lattice correlator eq. ()26|) in terms of these quantities. 
Let the point xq in the x-y plane denote the midpoint of the twisted plaquette Pxy{xo) 
with the weight CkU{pxy{xo)). The transfer matrix connects two such plaquettes in 
the z-direction (see fig. Q)). Then the correlator is: 




= Z-'Tt |yL.-r-iyi/2y^yWi/2 ^j.ik)y-^ Tf ^/^T^^Ty^} . (54) 

In this form we have an unwanted endpoint effect in the correlator, since the transfer 
matrices at the locations of the monopoles have only their inner potential factor 
twisted. The problem is remedied by using T' instead, in terms of which one simply 

gets 

q,)(r) = Z-'Tt {t'^-^ (t'(^'))'] . (55) 
Using eq. (j53|) we finally obtain our desired form for the correlator, 

C(fc)(r) = Z-^Tr {t'^^"^ V^T'^ V^} . (56) 

Since T and T' have the same spectrum, we thus conclude that the monopole pair 
correlator decays with distance, governed by the physical spectrum of the transfer 
matrix. For large r, the screening mass ttlm should thus represent the lightest state 
of the spectrum that couples to, and in Yang-Mills theory this would be the 0+"'" 
glueball. 

5.4 Generalization to 3+1 dimensions. 

A similar discussion applies to the case where a small periodic space dimension of 
length 1/T is added to our large two dimensional x — y space. The line of twisted 
plaquettes is then extended to a surface of twisted plaquettes which is closed in this 
periodic direction and defines the twisted action S(^k) in eq. (|TH|l . 

In the operator formalism our vortex operator is then repeated for every slice 
along a closed loop in this periodic direction and the product of all is denoted by 
Vk- Closure of this surface in the periodic direction guarantees that no plaquettes 
with a link in the periodic direction will be twisted. Choosing the transfer matrix T' 
of eq. (jH^ for the correlator of the extended vortex operators reproduces the path 
integral with the twisted action S'(fc). Finally, the length 1/T can be taken to infinity. 

In conclusion, in an infinite volume the correlator eq. ()26|l is measuring the lowest 
mass of the states excited by the vortex operator. This is true for any finite lattice 
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spacing. In SU(2) there is only one set of states excited, and the 0^ is the lowest one 
in the continuum limit. 



6 Lattice results in 3d 

The goal of our numerical study is to measure the screening mass tum from the 
large-distance behaviour of the monopole correlation Cfc(r), eq. (j2EI)- We work in the 
infinite temperature limit, where the system is effectively 3d and the temperature 
dependence of the free energy can be dropped, so 

Cik)ir) = = e-^-M/^, F„(r) = Fq - - exp(-mMr), (57) 

Z r 

(cf eq. (1) and (18)) corresponding to the exchange of a boson of mass rriM in 3 
dimensions. Z(fc)(r) is the partition function of a system with a stack of k twisted (or 
"flipped" for SU{2)) plaquettes as per eq.l^^. 

To compute ^^^^^ our strategy consists of factorizing the ratio into r factors, 
each of order 1, as done in Ref. [S] for the 4-dimensional case: 

Mr) = X Itnil X „ X £W. (58) 

Z Z{r-\) Z{r-2) Z ' ^ ' 

where Zir) = Z(fc)(r) and Z{j) has j twisted plaquettes. Each factor can be written 
as an expectation value: 

= (exp(|ReTr(C. - l)f/(j>,))^ , (59) 

where U{pj) is the j*'^ plaquette to be twisted, and the expectation value is taken 
with respect to Z{j — 1). Equivalently, the same ratio can be expressed as a ratio of 
two expectation values with respect to an interpolating partition function Z{j — 1/2), 
where U{pj) is multiphed by ^'y^: 



Zjj) (exp(|ReTr^[/(p,))),_i/2 



ZU - 1) (exp(|ReTr lf^^(p,))),-i/2' 



(60) 



This latter expression has a smaller variance than eq. ()59|) . Moreover, for SU(2) 
= 0, so that plaquette pj has zero coupling in Z(j — 1/2). Since the 4 links making 
up U{pj) are thus decoupled, one can form estimates of the two observables, from 
n estimates of each of the 4 links. This provides additional variance reduction. 
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On a lattice of size 32^ x 64, for each value of j = 1, 2, .., 32 (or in practice until 
the signal became undetectably small), we have performed independent Monte Carlo 
simulations, collecting 1 — 3 10^ measurements in each. Each simulation provides an 
estimate of z^ji^i^ = ciji-a) ^ different j. We can then compare all these estimates 
to the functional form 



log- 



C{r) 



Z]i=(f Cj/r exp(— m^r) 



C{r - a) X^i'o" Ci/(r - a) exp(-mi(r - a)) ' 
corresponding to the exchange of (zmax + 1) increasingly heavy scalar bosons of masses 



The successive ratios 



C(r) 



are shown in Fig. 01 A single boson exchange is 



C(r-a) 

insufficient to describe the data. Therefore, we fit the data to an Ansatz corresponding 
to the exchange of 2 bosons of masses mo and mi, over the range r G [3, 20]. 



100 



o 
o 




10 

x/a(P=16.7) 



Figure 3: Monte Carlo measurements of C{x)/C{x — a) (eq. ^U\i ) versus distance. 
Simulations at /3 = 16.7 and 12.0 give consistent results. A 2-mass fit per eq.(pT|) 
describes the data well, from x/a = 3 on. 

The fitted masses are shown in Table I. To check that we are simulating contin- 
uum physics, we repeated the simulations at two different values of the lattice spacing, 
corresponding to [3 = 16.7 and 12.0. The fitted masses from each simulation are con- 
sistent. Thus, we can fit both datasets together, obtaining the last line of Table I. The 
groundstate mass mo/gj^ = 1.69(8) is compatible with the mass of the 0"^^ glueball 
measured in the zero-temperature 3d SU(2) theory [12], m{0~^~^) / gl^ = 1.58(2). We 
interpret the fact that it is somewhat larger as being due to remaining contamina- 
tion of higher excitations (for example, no zero momentum projection is used in the 
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mo/gj4 




16.7 
12.0 


1.73(12) 
1.54(10) 


3.75(19) 
3.43(10) 


together 


1.69(8) 


3.67(12) 



Table 1: Mass of the hghtest exchanged boson (mo/gl^) and next-hghtest {mi/g^^) 
from simulations at two different values of the lattice spacing in the scaling region, 
corresponding to the curves Fig. |21 The two sets of results are consistent and can be 
fitted together (bottom line). 

monopole correlator.) On lattices of the sizes used here, glueball masses are known 
to be free of finite size effects, so we believe our infinite volume discussion is justified. 

7 Summary and prospects 

In this paper we have determined some simple but useful properties of the magnetic 
Z(N) screening mass. In particular, we have shown that the screening mass has 
the numerical value of the Oj;"*" state of a spatial Hamiltonian with a periodic space 
dimension of length 1/T. So the screening mass can be computed by correlating 
simple local operators with the right quantum numbers. For the extreme case of 
infinitely large T we have shown by lattice simulation that indeed the magnetic mass 
as measured by twisted actions and the 0^"^ state in 2+1 dimensions do coincide. 
The latter is very encouraging for the technique of simulating the twisted action in 
cases where one cannot revert to a simpler operator, as in ref. P|. 

It is very desirable to evaluate the magnetic mass for all T between the two cases 
where they are already known. How sensitive is the magnetic mass to the Z(N) 
transition that governs the behaviour of the thermal Polyakov line? One would say 
that the real part of the Polyakov line has the required quantum numbers to excite 
the 0^+. Thus one might be led to think that the magnetic mass should follow the 
behaviour of the Polyakov line: second order in SU(2), almost second order in SU(3), 
and first order for N larger than 3. However, Z(N) symmetry operators do not affect 
local states . Only below there are torelon states (i.e. string states winding in the 
periodic direction j2]), that do transform non-trivially, and indeed are disappearing 
above Tc- Ref. [H] shows that the screening mass smoothly increases with temperature, 
with no special feature at T^.. However, for SU(3) the mass of the representing the 
g++ ground state dips somewhat as one approaches Tc [12]. We thus feel the jury is 
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still out on the behaviour of magnetic screening and the spatial string tension near 
the critical point. 
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